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Abstract 

We examine the one loop contributions arising in the Two-Higgs-Doublet Model 
(THDM) to the W-boson anomalous magnetic dipole and electric quadrupole form 
factors for both photon and Z couplings relevant at collider energies. While the model 
parameter and (^-dependencies of these form factors are found to be significant, the 
corresponding size of these corrections are relatively small in comparison to unity. They 
are, however, found to be comparable in magnitude to the usual Standard Model loop 
corrections. Radiative corrections to the Higgs particle masses and couplings due to 
heavy top-quarks are included in the analysis. 



*Work supported by the U.S. Department of Energy, Division of High Energy Physics, Contracts 
W-31-109-ENG-38 and W-7405-Eng-82. 



Many theories which go beyond the Standard Model (SM), such as the Peccei-Quinn 
Model, Supersymmetry(SUSY), and E 6 sup erst ing-inspired models, require the scalar sector 
to be enlarged by the addition of an extra Higgs doublet (I], |3| . Such Two-Higgs-Doublet 
Models (THDM) have an extremely rich phenomenology 0] as the particle spectrum now 
contains two neutral CP-even states, h and H, a neutral CP-odd state, A, and a conjugate 
pair of charged fields, H 1 * 1 . Unfortunately, although there have been direct searches for 
such particles at both e + e~0 and hadron|J colliders as well as searches using indirect 
techniques^, they have failed to be observed. Of course, the searches so far conducted have 
only explored a small region of the allowed parameter space in these models but the advent 



of hadron supercolliders || |9| and e e machines in the TeV range |10| will allow essentially 
all of the parameter space of current interest to be probed. 

As the detailed predictions of the SM begin to be examined, one way to look for new 
physics is via loop-effects through possible small deviations away from SM expectations. 
Using existing data, this has already been done by a large number of authors [ITT]] for many 
possible SM extensions. As has been suggested [|12], one place to look for deviations in the 
future is the trilinear gauge couplings, 7WW and ZWW, which are uniquely predicted in 
the SM. In the 7WW case, with all particles on-shell, loop corrections to the tree-level SM 
predictions were considered some time agofTBf for several different SM extensions including 
the THDM case. In this paper, we return to a discussion of the loop corrections to the 7WW 
vertex within the context of the THDM, incorporating several refinements to fully cover the 
parameter space and make a direct connection to the physical Higgs spectrum. We then 
extend the previous analysis to the ZWW case and generalize to the situation where the 7 or 
Z can be off-shell. This is, of course, what will be realized experimentally when the trilinear 
couplings are probed in W-boson pair production at e + e~ colliders. In particular, we are 
interested in the g 2 -dependence of the anomalous magnetic dipole and electric quadrupole 
moment form factors and how they differ from the expectations of the SM. We find that for 
both the 7WW and ZWW cases, and all interesting values of q 2 , both the model parameter 
and g 2 -dependencies of these deviations are significant. Unfortunately, we also show that 
the size of these deviations are always relatively small in comparison to unity although they 
can be comparable in magnitude to the SM predictions themselves. It may, however, be 
possible to probe these anomalous couplings at the few parts per mil level at the NLCP^j 
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given sufficient integrated luminosity. 

To perform our calculation we label the momentum flow in the V(= 7 or Z)WW vertex 
as shown in Fig.l. Although we allow V to be off-shell, we assume that terms proportional 
to Q x can be dropped since V is assumed to couple to massless fermions. Both W's are 
assumed to be on-shell in the discussion below. Following the notation of Hagiwara et.al. 



12], the relevant part of the coupling in momentum space is given by 



v 



Mi 



w 



when both the W + and W are on-shell. At tree-level, the SM gauge symmetries demand 

that fY = 1, fX = 0, and fY = 2 with g^ww = — e and gzww = ~9Cw where c w = cos9 w . 
At one-loop, however, these parameters are now shifted and are given instead by (here s 
is the square of V invariant mass, i.e., s = q 2 , which we will assume, for the moment, is 
positive) 

fY = Ay (2) 
fY = 2 + AgY + AK V + X v 

so that the momentum space coupling can be written as 



ff^rjr = (i + Asr + ^A y )[2A^ + 4(QV A - QTsr x )] (3) 

- 8^rP X Q»Q u + 2[Ak v + X v (l - s/M*y)](Q"g>* - Q»g vX ) 

M w 

Note that both Ay and Any are finite and g 2 -dependent, whereas AgY contains ultra-violet 

infinities corresponding to a charge renormalization. (In the case where a massless scalar 
is exchanged between the W's, an infra-red divergence can also be isolated in AgY which 
is then cancelled by a real emmission diagram.) Instead of being constants, both Any and 
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Ay are true form-factors. As can be seen from this equation, the form of the finite terms 
that we are interested in here (Any and Ay) are somewhat different than those used by 
Couture et.al.|l3| in the case where the initial photon, as well as the W-pair, are on-shell. 



This problem is easily overcome by noting that the two parameterizations are simply related 
via 

(AK) C outure = Ak V + Ay(l - s/M^) 

(AQ)couture = ~ 2Ay (4) 



Our convention is now the standard one adopted by most authors in the literature [12 . 



The two classes of diagrams contributing to Any and Ay at one loop in the THDM 
are shown in Fig. 2. In the first class, the external V couples to either a pair of W's, a pair 
of charged Goldstone bosons, or one of each and the two external W's are connected by the 
exchange of either h or H. In diagrams of the second class, the external V couples to H 1 * 1 and 
the external W's are connected by h, H, or A exchange. (Other classes of diagrams involving 
the ZhA and ZHA vertices for the V=Z case are found to cancel due to Bose statistics.) 



Denoting the contributions of these twos set of diagrams by the q -dependent quantities I\ 
and I2, we can write the entire THDM contribution to either Any or Ay symbolically as 

\JI x {h) + I 2 (H, F±)]4_ a + ifh(H) + I 2 (h, F±)]4_ Q + I 2 (A tf±) - SM (5) 

In this equation, the aguments of l\ i represent the relevant particle masses, cp_ a = cos(f3—a) 

and S/3_ a = sin{(3 — a) which are the usual mixing angle factors^], and / depends on the 
choice made for V: 

f 1, ^ = 7 

[ {\-2x w )/2(\-x w ), V = Z 

Also note that in Eq.(5) we have subtracted off the SM contribution to both Afty and Ay 

since we are interested in how much these quantities are changed when the SM is extended 
to the THDM case. To be specific, we have subtracted the SM contributions assuming a SM 
Higgs mass of 200 GeV and taking the appropriate value of q 2 . (Of course, a different choice 
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for the SM Higgs mass will only shift our results for Any and Ay by a pair of constants.) 
This subtraction will be most noticable in the case of Ay which would normally vanish in the 
limit that all the masses in the loop become large. Instead, we simply recover the negative 
of the SM contribution in the present case. Note that in the limit that H, A, and H 1 * 1 masses 
get very large and s@- a — * 1, we will recover the SM limit provided that h is identified with 
the Higgs scalar of the SM in which case the difference defined in Eq.(5) will vanish. 

Since the THDM contains many unknown quantities, it is very difficult to explore all 
of the parameter space. The above results, while being valid in any THDM whose scalar 
sector conserves CP, are difficult to analyse numerically. In order to simplify our results and 
to reduce the number of free parameters we have assumed that the THDM Higgs masses 
and couplings are as given by SUSY- version of the Model with the one-loop heavy top-quark 
mass corrections included fl5|]. Specifically, we have chosen the same input parameters as 
used in the effective potential analysis of Barger et.al. || except that we have taken the 
top-quark mass to be 130 GeV. Our results are not qualitatively sensitive to this particular 
choice of parameters. With this great simplification, the only remaining free parameters are 
the ratio of the two Higgs doublet vacuum expecation values, tan/3 = v t /vb, and the mass of 
the CP-odd field A, itla- All other masses and mixing angle parameters are then fixed. For 



numerical purposes, will we assume My/ = 80.14 GeV |16| and Mz = 91.175 GeV |17| with 
cw = M w /M z in the analysis that follows. The functions Ii )2 can be expressed generically 
in terms of a small set of integrals which reduce to those of Couture et al. ||13|| in the limit 



that q 2 — > 0. For example, in the Ay case, we define the integral 

J(r q ,r p ,r ) = 6a f 1 dt t 3 (l-t) f 1 du u(l - u)\A - Bu{l - u)}' 1 (7) 
Jo Jo 

with B = r q t 2 and A = t 2 — t(l + r — r p ) + r . The r's are given by 

r q = q 2 lM 2 w 

r p = ml ± /M 2 w (8) 
r = ml/M 2 , 

with m being a generic mass for a neutral Higgs field. Note that we have defined the 
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dimensionless parameter, a, as 

~ 96tT 2 

which sets the typically small scale for the size of the loop-induced anomalous couplings 
now find that Ji j2 in Eq.(5) can now be simply expressed as 

I?(h,H) = J(r q ,l,r h>H ) 
iSih^A-H*) = J(r q ,r p ,r KHA ) (10) 

in obvious notation. A similar set of relations can be developed for the case of Ak v ; definining 
the integrals 

K(r q ,r p ,r ) = —6a dt du t(3ut — 1) log\A — Bu(l — u)\ 

Jo Jo 

L(r q ,r p ,r ) = 3a jf* dt J* du t 2 (A + 2tr - 2r )[A - Bu(l - u)]' 1 (11) 
with A and B as given above, we find that 

I? K (h,H) = K(r q ,l,r h>H )-(l-r q )J(r q ,l,r h , H ) 
l2 K (h, H, A; i? ± ) = K(r q ,r p ,r h , HA ) - (l-r q )J(r q ,r p ,r hiHjA ) + L(r q ,r p ,r h , HA ) (12) 

The above integrals are relatively easy to perform numerically once care is given to potential 

imaginary parts. Since Ak v and Ay are themselves generated only at the one-loop level, 
they contain no imaginary parts at this order; such imaginary parts will, however, arise at 
the two-loop level. (The imaginary part of the one-loop vertex function is in fact buried 
in the quantity AgY introduced above.) These results are seen to simply reduce to those 
obtained earlier assuming V=7 and taking q 2 to zero. 

In order to get a feeling for the relative size of the THDM contributions we note that 
the SM with m t = 130 GeV and M H = 200 GeV predicts A/t 7 = 13.5a and A 7 = -0.6a 
for q 2 = 0. (Note that the predicted size of A 7 is more than an order of magnitude smaller 
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(9) 
. We 



than that of Ak 7 .) As we will see, the THDM contributions are quite comparable to these 
values in size. Unfortunately, since a is a small quantity, A/t 7j z will almost always be less 
than ~ 0.01 in magnitude (with A 7i z even smaller) in the THDM case. It is important 
to remember that such SM contributions are subtracted out in what follows as purely SM 
radiative corrections are already well understood in processes such as e + e~ — > W + W~ in 
which our anomalous couplings might be measured. 

Since the q 2 = situation has already been discussed previously for the case V — 7 
H| we will mention it here only to illustrate the q 2 dependence of the two form factors. 



It will, however, show many features in common with cases where q 2 7^ and where V=Z. 
Figs. 3a and 3b show Ak 7 and A 7 as functions of tua for different choices of tanf3 assuming 
q 2 =0. The first thing we observe is the apparent lack of sensitivity to the choice of tan/3, 
a property that will persist in all other cases as we will see below. Both Ak 7 and A 7 have 
their greatest magnitudes when tua is small. The reason for this is clear; as itla increases 
the corresponding masses for H and H 1 * 1 also become large leaving only h light. The h 
contribution is then largely cancelled by the usual SM subtraction; this is particularly true 
in the A 7 case where decoupling in the heavy mass limit is most effective. We also note 
again the general feature that the relative sizes of A/c 7 and A 7 generally differ by more than 
an order of magnitude; this effect will occur in all the cases we will examine below. As 
advertised, the THDM contributions to both of these parameters are generally comparable 
to their SM values in magnitude, particularly for small values of w,a- These results were 
found to be only moderately sensitive to the existence of heavy top-quark loop corrections 
to the Higgs' masses and couplings, particularly, at smaller values of tua, and only midly 
sensitive to the particular SUSY parameter values used to calculate these corrections. These 
properties continue to valid in all of the results we present below. 

To verify the lack of tanfi sensitivity, we show in Figs. 3c and 3d the values of Ak 7 and 
A 7 as functions of tanfi for different values of tua assuming q 2 =0. Quite generally, we see 
that the values of these parameters are approximately constant for tanj3 larger than unity 
which is just the region in which we expect the SUSY THDM mass and coupling realtions 
from the effective potential approach that we have used in these calculations to be most valid 
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What happens when we increase q 2 and we are no longer probing the static moments? 
Figs.4a-d and 5a-d show the results akin to Figs.3a-b for different values of q 2 and we see a 
significant change in the shape of the curves when the W-pair threshold is exceeded. The 
general lack of sensitivity to tan(3 is still evident except for some very small ranges of m^; 
this we have examined explicitly for several fixed values of q 2 . In the case of A 7 , we note 
that the magnitude decreases significantly as q 2 increases. This behaviour might have been 
anticipated by the results in Eq.(lO). Ak 7 also decreases in magnitude once the value of q 2 
exceeds the W-pair threshold but does not decrease significantly beyond v /s=200 GeV. To 
explore this behaviour further, we display the explicit g 2 -dependencies of both Ak 7 and A 7 
in Figs. 6a-b assuming tan(3 =2 for different values of tua- (Due to the general overall lack 
of sensitivity to tanf3, these results are largely independent of the value of this parameter.) 
Except for the region near the W-pair threshold, these functions are generally smooth and 
A 7 is seen to tend to zero in the large q 2 limit as expected. The asymptotic behavior of Ak 7 
is somewhat more complex in this limit. We do, however, see the general result that the 
(^-dependence of both Ak 7 and A 7 in the THDM cannot be ignored in that the deviations 
expected for the static moments in the THDM are quite different from those one would 
measure at a high energy collider producing W-pairs. Unfortunately, although the size of 
these corrections we have found are comparable to the SM results themselves, they are 
numerically quite small. 

Do these general features persist when we look at the case V=Z as we might nat- 
urally expect? Fig. 7 answers this question in the affirmative; here we see repeated lack of 
dependence of our results on tanf3 and the drastic reduction in magnitude of Xz as 
increases. Generally, the magnitudes we find for Xz and Akz are comparable to what was 
seen in the case V=7. In fact, the numerical results for the two quantities in the V=Z and 
V=7 cases are found to be very similar for most values of the parameters. Fig. 8 shows that 
the g 2 -dependence of the form factors in both cases is also extremely similar as one might 
expect. Thus in the case of W-pair production at ^=200 GeV, the anticipated size of A 7 
amd Xz (as well as Ak 7 and Ak z ) are quite comparable in both the SM and the THDM. 

We have examined the contributions to the anomalous WWZ and WW7 trilinear 
couplings within the context of the THDM over a wide range of model parameters and 
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values of q 2 . We have incorporated the one- loop radiative corrections due to a heavy top- 
quark into the mass and coupling relationships amongst the various physical Higgs fields. 
The major results that we have obtained in our analysis are as follows: 

(i) Both Any and Ay were found to be relatively insensitive to the choice of tanf3 
once this quantity was greater than unity The sensitivity to the value of m A was particularly 
stong especially for m A values below about 200 GeV. The usual SM results were obtained in 
the limit of large m A . Although the THDM contributions did not lead to order of mamnitude 
changes in the anomalous moment for factors, we did see that the additional contributions 
could be quite comparable to the SM results. 

(ii) The anomalous couplings were seen to be of comparable magnitude in both the 
V=Z and V=7 cases. Thus it would be impossible to seperate these two sources of anomalous 
couplings using e + e~ data alone. Although the values we have obtained for the anomalous 
moment form factors in the THDM remain reasonably small, measurements taken at the 
NLC may help to constrain the THDM parameter space. 

(Hi) The two form factors we examined were observed to exhibit a strong q 2 depen- 
dence. Their values, as would be measured via W-pair production at an e + e~ collider, are 
substantially different than their static values obtained at q 2 =0. \y was found to decrease 
in magnitude quite rapidly with increasing q 2 . 

(iv) The influence of incorporating the heavy top-quark loop corrections to the Higgs' 
masses and couplings was found to be significant particularly in the range of small values. 
The specific results were, however, quite insensitive to the particular choice of the SUSY 
parameters that were employed. 

Precision measurements of the anomalous moment couplings of the W at e + e~ col- 
liders and elsewhere may lead to the first indications of the existence of new physics beyond 
the Standard Model. 
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Figure Captions 



Figure 1. Feynman diagram for the VWW vertex with the momenta labeled as used in the text. 

Figure 2. The two classes of triangle diagrams responsible for generating the anomalous moments 
of the W-boson in the THDM. 

Figure 3. The 7WW anomalous moments for the case q 2 =0. (a) Ak 7 as a function of m A for 
tanfl =1, 2, 5, 10, and 30; (b) same as (a) but for A 7 ; (c) same as (a) but as a function of 
tanf3 for m^=50 (top solid), lOO(dashed-dot), 250(dashes), 500( dots), or 1000(bottom 
solid) GeV; (d) same as (c) but for A 7 with the top solid curve now corresponding to 
mA=1000 GeV and the order of the subsequent curves in (c) now reversed. 

Figure 4. m A dependence of A/t 7 for y/s= (a) 100; (b)200; (c)500; or (d) 1000 GeV with tan/3 as 
given in Fig. 3. 



Figure 5. Same as Fig. 4 but for A 



7" 



Figure 6. y/s dependence of (a) A/t 7 and (b) A 7 assuming tan(3=2. In (a), from top to bottom 
on the left hand side, the curves correspond to m A =50, 100, 200, 500, or 1000 GeV 
while in (b) the same order is reversed. 

Figure 7. ttia dependence of Akz and \z for the same values of tanf3 as in Fig.3 with either y/s 
=200, 500, or 1000 GeV. 

Figure 8. Same as Fig. 6 but for Ak z and \ Z - 
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